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Abstract. We present a detailed empirical assessment of how the galaxy luminosity function and stellar luminosity den- 
sity evolves over the last half of the universe's age (0.2 < z < 1.2) for galaxies of different spectral energy distributions 
(SED). The results are based on ~ 25000 galaxies (R ^ 24) with redshift measurements (ctz ~ 0.03) and SEDs across 
Aobs ~ 350 . . . 930 nm. The redshifts and SEDs were derived from medium-band photometry in 17 filters, observed as part 
of the COMBO-17 survey ("Classifying Objects by Medium-Band Observations in 17 Filters") over three disjoint fields with a 
total area of 0.78 square degrees. Luminosity functions (LF), binned in redshift and SED-type, are presented in the restframe 
passbands of the SDSS r-band, the Johnson B-band and a synthetic UV continuum band at 280 nm. 

We find that the luminosity function depends strongly on SED-type at all redshifts covered. The shape of the LF, i.e. the 
faint-end power-law slope, does depend on SED type, but not on redshift. However, the redshift evolution of the characteristic 
luminosity A/* and density (j)* depends strongly on SED-type: (1) Early-type galaxies, defined as redder than a present-day 
reference Sa spectrum, become drastically more abundant towards low redshift, by a factor of 10 in the number density (p* from 
z=l.l to now, and by a factor of 4 in their contribution to the co-moving r-band luminosity density, j^. (2) Galaxies resembling 
present-day Sa- to Sbc-colours show a co-moving number density and contribution to jr that does not vary much with redshift. 
(3) Galaxies with blue spectra reflecting strong star formation decrease towards low redshift both in luminosity and density, 
and by a factor of 4 in their jV contribution. Summed over all SED types and galaxy luminosities, the comoving luminosity 
density decreases towards low redshift, between z=l.l and now, by a small amount in restframe r and B, but by a factor of ~ 6 
in restframe 280 nm. Al z = 1.1, galaxies redder than She's, contribute 40% to the total jr, which increases to 75% by z=0. 
For Arost =280 nm, this increase is from 12% to 25% over the same redshift interval. 

Comparison of the three independent sight-lines shows that our results are not significantly affected by large-scale structure. 
Our lowest redshift bin at z — [0.2, 0.4] largely agrees with the recent assessment of the present-day galaxy population by 
SDSS and 2dFGRS and deviates only by an excess of "faint blue galaxies" at 2 ~ 0.3 compared to very local samples. Overall 
our findings provide a set of new and much more precise constraints to model the waning of overall star formation activity, the 
demise of star-bursts and the strong emergence of "old" galaxies, with hardly any young population, over the last 6-8 Gigayears. 
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1. Introduction 

The formation and subsequent evolution of galaxies are deter- 
mined both by the overall gravitational growth of structure and 
by the physics of gas cooling, star formation and feed-back 
which determine the successive conversion of gas into stars. 
Hierarchical structure formation within cold dark matter sce- 
narios and their various extension to address star formation, 
now provide a comprehensive, but parameterized, framework 
for galaxy formation (e.g. Cole et al. 2000). 

The onset of galaxy formation seems to take place at such 
high redshifts that it has so far escaped direct observation. 



Instead, most observations of galaxies have so far concentrated 
on obtaining large samples at lower, more easily accessible red- 
shifts. 

At the high redshift end, pioneering galaxy survey work has 
meanwhile reached redshifts z ^ 5 (e.g. Ouchi etal. 2002). In 
the local universe, two large surveys are currently characteriz- 
ing in detail the luminosity function of the galaxy population, 
based on large samples with > 10^ objects, the 2dF Galaxy 
Redshift Survey (2dFGRS, e.g. Madgwick et al. 2002) and the 
Sloan Digital Sky Survey (SDSS, e.g. Blanton et al. 2001). 

During the past ten years several studies have aimed to map 
out the evolution of the luminosity function and the total lu- 
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minosity density from the local universe to a redshift of ~ 1 



(Lilly et al. 1995 




Madau, Pozzetti & Dickinson 1998 


; Lin 


etal. 1999; 


Fried etal. 2001 


). But samples sizes well in ex- 



cess of a few thousand objects are only becoming available now 
or in the near future, e.g. with the 17-colour survey COMBO- 
17 presented here and with the large spectroscopic campaigns 
DEEP (iKoo 20011 Im et al. 20021) and VIRMOS (LeFevre 



2001). The scientific inferences from existing, faint surveys out 
to z ^ 1 have been limited mainly by their sample sizes, ag- 
gravated by the strong influence of large-scale structure when 
observing small co-moving volumes. The present survey, and 
other ongoing initiatives, aim at improving the measurement 
of the luminosity function by smoothing over structure and in- 
creasing the volume. 

The COMBO-17 project ("Classifying Objects by 
Medium-Band Observations in 17 Filters") was designed to 
provide a sample of ^^50,000 galaxies and < 1,000 quasars 
with rather precise photometric redshifts based on 17 colours. 
In practice, such a filter set provides a redshift accuracy of 
^2, gal « 0.03, cTz.QSO ^ 0.1, smoothing the true redshift 
distribution of the sample only slightly and allowing the 
derivation of luminosity functions. 

The foremost data analysis goal of the COMBO-17 ap- 
proach is to convert the photometric observations into a very- 
low-resolution spectrum that allows simultaneously a reliable 
spectral classification of stars, galaxies of different types and 
QSOs as well as an accurate redshift (or SED) estimation for 
the latter two. The full survey catalogue should contain about 
75,000 objects with classifications and redshifts on 1.5 of 
area. This fuzzy spectroscopy consciously compromises on red- 
shift accuracy (a^ w 0.03) in order to obtain very large sam- 
ples of galaxies with a reasonable observational effort. While 
both characteristics are well suited for the analysis of an evolv- 
ing population, they understandably do not permit dynamical 
or chemical studies which require quite detailed spectroscopic 
information. 

While the photometric redshift technique has akeady been 



applied to galaxy samples about 40 years ago ( |Baum 1963 
Butchins 1983), we have optimized the technique by in- 



creasing the number of filters and narrowing their bandwidth 
to obtain better spectral resolution and more spectral bins. 
Therefore, COMBO-17 also provides identifications and rea- 
sonably accurate redshifts for quasars (see also Koo 1999 for 
a nice overview on photometric redshifts and SubbaRao et al. 
1996 on applying the technique in the context of luminosity 
functions). 

The goal of the present paper is the use of COMBO-17 red- 
shifts and SEDs for 25,000 galaxies over 0.78 to draw up 
a detailed, empirical picture of how the population of galaxies 
evolved over the last half of tiie universe's age. 

Our paper is organized as follows: in Sect. 2 we present the 
observations that have led to the current sample of ^ 25,000 
galaxies. Our techniques for obtaining their redshifts, SED 
classification, luminosities and completeness are described 
in Sect. 3. The resulting sample properties are discussed in 
Sect. 4. In Sect. 5 we derive our "quasi-local" luminosity func- 
tion, drawn from the redshift interval of z = [0.2, 0.4] and com- 
pare it with the results from more local samples obtained by the 



2dF Galaxy Redshift Survey and the Sloan Digital Sky Survey. 
Finally, we show the evolution of the luminosity function and 
the luminosity density out to z < 1.2. 



2. Observations - the COMBO-17 survey 

The COMBO-17 survey has produced multi-colour data in 17 
optical filters on 1 of sky at high galactic latitudes, including 
to date the Chandra Deep Field South (CDFS) and the field of 
the supercluster Abell 901/902. The filter set (Fig. |l]and Tab. |l|) 
contains five broad-band filters (UBVRI) and 12 medium-band 
filters stretching from 400 to 930 nm in wavelength coverage. 

All observations presented here were obtained with the 
Wide Field Imager (WFl, Baade et al., 1998, 1999) at the 
MPG/ESO 2.2-m telescope on La Silla, Chile. They encom- 
pass a total exposure time of ^^160 ksec per field including a 
^^20 ksec exposure in the R-band with seeing below 0'.'8. The 
WFl provides a field of view of 34' x 33' on a CCD mosaic 
consisting of eight 2k x 4k CCDs with ~ 67 million pixels 
providing a scale of O'.'238/pixel. The observations started in 
the commissioning phase of the WFl in January 1999 and are 
continuing as the area is extended to cover more fields. 

The instrument design and the survey concept have been 
matched to the requirements of deep extragalactic surveys. The 
morphological and spectral dataset of COMBO-17 is primarily 
intended for studies of (i) gravitational lensing and (ii) evo- 
lution of galaxies and QSOs. Indeed, the optics of the instru- 
ment have been designed with the prime application of lensing 
in mind, while the filter set was tailored for the task of object 
classification and redshift estimation. 

Observations and data analysis have been completed for 
three fields (see Table ^ covering an area of 0.78 D° and pro- 
viding a catalogue of ^200,000 objects found by SExtractor 
( Bertin & Arnouts 1996[ ) on deep, high-resolution R-band im- 
ages with 5(7 point source limits of i? « 26. 

These deep R-band images provide very sensitive surface 
brightness limits and allow to establish the total object photom- 
etry using the SExtractor measurement MAG-AUTO. Except 
for L-stars and quasars at z > 5, they provide the highest 
signal-to-noise ratio for object detection and position measure- 
ment among all data available in the survey. 

The spectral shapes of the objects in the R-band se- 
lected catalogue were measured with a different approach. 
Photometry was obtained in 17 different passbands by project- 
ing the object coordinates into the frames of reference of each 
single exposure and measuring the object fluxes at the given lo- 
cations. In order to optimize the signal-to-noise ratio, we mea- 
sure the spectral shape in the high surface brightness regions of 
the objects and ignore potential low surface brightness features 
at large distance from the center. 

Since seeing variations among the different bands would 
introduce artificial colour offsets by changing observing con- 
ditions typical for ground-based observations, we need a non- 
standard photometry approach to measure spectral shapes ac- 
curately. In fact, we need to measure the same central fraction 
of an object in every band as it would appear in equal seeing. 
To this end, we employ a seeing-adaptive, weighted aperture 
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Table 1. The COMBO-17 filter set: Exposure times and lOcr 
magnitude Hmits reached for point sources, averaged over all 
three fields. The i?-band observations were selected to be taken 
under the best seeing conditions (FWHM 0'.'55 . . . O'.'S). 
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photometry as performed by the package MPIAPHOT (Roser 



& Meisenheimer 1991; Meisenheimer, in prep.) 



MPIAPHOT measures the central surface brightness of ob- 
jects after convolving their appearance outside the atmosphere 
to an effective PSF of I'.'S diameter. In detail, the procedure 
measures the observed stellar PSF on each individual frame 
and chooses the necessary Gaussian smoothing for reaching a 
common effective PSF of I'.'S uniformly on all frames in all 
bands. For most objects this measurement is similar to a flux 
measurement in an aperture of ~ 2" diameter in 1'.'5 seeing. 

The photometric calibration is based on a system of faint 
standard stars in the COMBO-17 fields, which we established 
by spectrophotometric calibration with respect to spectrophoto- 
metric standard stars in photometric nights. Our standards were 
selected from the Hamburg/ESO survey database (Wisotzki et 
al. 2000) of digital objective prism spectra (see Wolf et al., 
2001b, for procedure). By having standard stars within each 
survey exposure, we were independent from photometric con- 
ditions for imaging. 

In summary, all luminosities used in the paper are based 
on SExtractor MAG- AUTO measurements on the deep R-band 
stack, while all redshift and SED fits are based on seeing- 
adjusted aperture measurements (with MPIAPHOT) across all 
bandpasses. For more details of the data reduction, we like to 
refer the reader to a forthcoming technical survey paper (Wolf 
et al. in prep.). In this paper, all magnitudes are cited with ref- 
erence to Vega as a zero point. 

3. The galaxy catalogue 

The galaxy catalogue is extracted from the full survey cata- 
logue purely on the basis of spectral information. There are 
no morphological criteria used to differentiate between stars, 
galaxies and quasars. Indeed, many faint galaxies appear com- 



Table 2. Positions and galactic reddening (Schlegel et al. 1998) 
for the three COMBO-17 fields analysed. All observations 
were obtained at the Wide Field Imager at the MPG/ESO 
2.2m-telescope at La Silla. 



Field 



«J2000 



gal 



CDFS 
A 901 
S 11 



03'' 32" 



25' 
17" 
58" 



-27° 48' 50" 
-10°01'25" 
-01°42'50" 



223? 6 
248? 
270? 5 



-54?5 0.01 
+33? 6 0.06 
+56? 8 0.02 



pact in typical ground-based seeing, while binary stars can 
produce objects with stellar spectra but extended appearance. 
Therefore, the abundance of photometric information provides 
a safer separation between the object classes if analysed with a 
classification technique as presented in the following section. 

3. 1. Classification and redsliift estimation 



The photometric measurements from 17 filters provide a low- 
resolution spectrum for each object to be analysed by a statisti- 
cal technique for classi fication and redshift estimation base d 
on spectral templates ( Wolf, Meisenheimer & Roser 2001 ). 
This approach has already been applied to the Calar Alto Deep 
Imaging Survey (CADIS) and to a mass determination of the 
galaxy cluster Abell 1689 ( [Wolfet al. 20014 |Dye et al. 2000| ). 

Since these initial analyses, we have improved the galaxy 
templates in the restframe UV region where the original Kinney 
templates E, SO, Sa and Sb show fairly noisy patches in the 
wavelength interval of A = 290 . . . 330 nm. These patches have 
been replaced by spectra obtained with the stellar population 



synthesis code PEGASE (Fioc & Rocca-Volmerange 1997) by 



first matching these to the Kinney templates. For the purpose 
of efficiency, we have also changed the redshift axis in the grid 
by making it equidistant on a log{l + z) axis rather than on 
a linear z axis. We have not yet incorporated trustworthy tem- 
plate information bluewards of the Lyman-alpha line, and still 
restrict the redshift range such that the existing templates al- 
ways cover our entire filter set. This constrains the investigated 
redshifts to z < 1.55 for now and leads to a deliberate exclu- 
sion of higher redshift objects from the catalogue. While they 
are not within the scope of this paper, some high-redshift galax- 
ies, e.g. Lyman-break objects at z ~ 3 ... 4, could be mistaken 
by the redshift estimation to reside at low-redshift and contam- 
inate the sample to a very small degree at the faintest end. 

The quasar library has also been improved by deriv- 
ing it from the more modern SDSS QSO template spectrum 
(Budavari et al. 2001 ) rather than from the Francis et al. (1991) 



emission line contour. This will make a difference mainly for 
a better detection of low-redshift quasars and hardly affect the 
rich class of galaxies in a statistical sense. Likewise the stel- 
lar library has been improved by omitting stars of spectral type 
O, B and A from the Pickles (1998) atlas, which we do not 
expect to see in the Galactic halo anyway, while template spec- 
tra of white dwarfs, subdwarfs and blue horizontal branch stars 
have been included. For all details we like to refer the reader to 
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Fig. 1. COMBO-17 filter set: Total system efficiencies are shown in the COMBO-17 passbands, including two telecope mirrors, 
WFI instrument, CCD detector and average La Silla atmosphere. Combining all observations provides a low-resolution spectrum 
for all objects in the field. Photometric calibrations of such "multi"-colour datasets are best achieved with spectrophotometric 
standards inside the target fields. 



our forthcoming pa per on the accuracy of the classification and 
redshift estimation ( Wolfetal. in prep ). 

The galaxy catalogue of concern in this paper is based on 
the observed average templates from Kinney et al. (1996), ex- 
cept for the modifications mentioned above. These ten tem- 
plates cover typical local galaxy SEDs from elliptical galaxies 
to starbursts (see Wolf, Meisenheimer & Roser, 2001, for de- 
tails). Altogether, they probably encompass the widest range 
of average ages possible for stellar populations in galaxies, 
but they explicitly do not contain templates for deeply dust- 
enshrouded starbursts as they may be part of the ERO galaxy 
population believed to reside at z — 1 ... 2. If galaxies with 
restframe SEDs similar to that of dust-enshrouded EROs were 
contained in this catalogue, they would not be identified as 
such, but rather with the best-fitting SED type among the ten 
Kinney templates, i.e. as likely old populations typical of lo- 
cal ellipticals. However, from an analysis of the CADIS galaxy 
sample we know that EROs are no t sufficiently abundant to 
change the conclusions of our study ( Thompson et al . 1999 ). 

Our classification and redshift estimation provides full 
probabilities for every object class based on all photometric 
measurements, rather than merely searching for a single best- 
fitting template. We therefore assign to an object the class 
which yields the largest sum over all probabilities for each 
spectral fit within that class. This prevents unreliable assign- 
ments of any class which may have a single spurious highly 
probable template, but otherwise fits very poorly over all other 
templates. While such a simple -minimisation is a powerful 
technique when plenty of spectral information is available (e.g. 
500 channels), the analysis of low-dimensional colour space 
benefits significantly from a wholesome probability-based ap- 
proach which handles ambiguities better (which is even more 
relevant in broad-band surveys with only five filters). 

While discriminative power is an important concern in ev- 
ery classification problem, completeness is no less important 
to avoid missing relevant fractions of a population. Therefore, 
it is reassuring, that less than 1% of the objects observed in 
COMBO-17 appear to have peculiar spectra which do not re- 
semble any template, but instead are outliers at more than a 



Table 3. Galaxy spectral types, their corresponding template 
range in the sequence of Kinney et al. (1996), and COMBO-17 
sample sizes. 



TYPE 


Template range 


SED type 


objects at 


objects at 




(Kinney et al.) 


range 


z =0.2-0.4 


z =0.2-1.2 


1 


E-Sa 


0-30 


344 


1365 


2 


Sa-Sbc 


30-55 


986 


5489 


3 


Sbc-SB6 


55-75 


1398 


7520 


4 


SB6-SB1 


75-99 


2946 


11057 


all 


E-SBl 


0-99 


5674 


25431 



3 (T level of significance. Among these are interesting individ- 
ual objects, but also a few blended objects and stars with un- 
corrected short-term variability skewing the observed spectra 
beyond our control. 

3.2. Galaxy spectral types 

The just described photometric data analysis of COMBO-17 
provides luminosity and restframe spectral type as the main ob- 
servables of galaxies. As a consequence, the evolution of their 
luminosity function can be investigated in sub-samples split by 
restframe spectral types. Note that, in contrast, types of mor- 
phology or explicit star formation history are not subject of the 
discussion presented here. 

There is no unique definition of SED types along a spectral 
parameter axis, since galaxies cover a continuum of parameter 
values. This fact not only applies to our own parameter defini- 
tion which is based on the grid of templates from Kinney et al. 
( 1 996), but is found in every definition of spectral type, whether 
it is based on restframe continuum colour as in the SDSS anal- 
ysis or derived from a principal component analysis (PCA) of 
the spectra themselves as in the 2dFGRS, an equivalent width 
of Ha-emission, or a mean age of a stellar population. We il- 
lustrate our choices and definitions in Fig. ||and Tab. |3] 

Having obtained very-low-resolution spectra from 17 pass- 
bands in COMBO- 1 7 it is only natural to use all available infor- 
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Fig. 2. Restframe colour of spectral templates and type defini- 
tion: 

Left-hand side: Restframe colours calculated for the sequence 
of galaxy templates along the entire SED parameter range. 
Shown are the restframe colours (m28o — B) vs. [B — r). The 
ten original Kinney templates are large labelled dots. The tem- 
plate sequence was mapped onto a grid of 100 templates by 
linear interpolation in flux space. 

Right-hand side: Restframe colours measured for the galaxy 
sample at i? < 23 and z = [0.2, 0.4] as obtained from the 



17-filter spectra (see Sect. 3.3). The four different SED types 



used in this paper are rendered in different colours. Since the 
SED types are determined from 17 passbands, they do not cor- 
respond to precisely defined intervals of any single colour axis. 



mation for the determination of spectral types, which is similar 
to the approach taken by the 2dFGRS. In fact, our classification 
procedure is based on a single-parameter set of templates and 
provides SED parameter estimates just as it provides redshift 
estimates. 

The range of SED parameter values just maps the Kinney 
et al. (1996) templates onto a sequence of values from (E 
galaxy template) to 99 (SBl starburst template), as in Wolf, 
Meisenheimer & Roser (2001). The distribution of SED values 
in the sample shows no obvious structures suggesting particular 
bins (see Fig. For the purpose of this paper, we decided 
to use a set of four SED types, derived by eye inspection of 
the evolutionary patterns observed, aiming at producing a clear 
picture of differential behaviour between the types (see Table ^ 
for the relation). Fig. || shows the mean templates representing 
the four types of galaxies at 0.2 < z < 1.2. 

The relationship between these types and the restframe 
colour indices formed among the passpands in our luminos- 
ity definition can be seen in Fig. ^. It shows that our template 
sequence (left panel) covers the colour distribution of the ob- 
served galaxies (right panel). It also illustrates that the SED 
classification (based on 17 filters) reflects a single-parameter 
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Fig. 3. Average templates for galaxy types: The average rest- 
frame templates of the four SED types at 0.2 < z < 1.2 are 
plotted with an offset of one flux unit each. The average tem- 
plate in typel becomes slightly redder from 2 ~ 1.1 to z ^ 0.3, 
but we see no change in the other types. The high and low red- 
shift versions of the type 1 template are plotted in grey, and the 
mean over all redshifts is shown in black. 



sequence, but does not correspond to precisely defined limits 
on a single colour axis, due to (a) the presence of noise, (b) 
the not entirely monotonic behaviour of the template sequence, 
and (c) the fact that the SED parameters of individual objects 
are determined from all 17 passbands and not from a single 
colour index. 

In comparison, the 2dFGRS analysis by Madgwick et al. 
(2002) measures a spectral type from a principal component 
analysis of all present spectra. An individual galaxy is then as- 
signed an 77 value, corresponding to the main parameter arising 
from the PCA. Much hke in COMBO-17, information from the 
whole spectrum was mapped onto a single parameter defining a 
sequence of spectral types, which could presumably be mapped 
onto our SED axis. The sequence was then split into four bins 
with the first one covering a prominent peak of no-emission- 
line galaxies in the 77-histogram. However, we did not attempt 
to reproduce their types since we lack the necessary colour data 
to perform such a task. 

On the contrary, the type definition in the SDSS analysis by 
Blanton et al. (2001) was based on a single colour index. The 
spread in reconstructed restframe (g — r) -colour was split into 
five bins of equal width on the magnitude scale. This is roughly 
equivalent to drawing equidistant vertical lines into Fig. || for 
separating the subsamples along their restframe (i3 — r) -colour. 
The figure demonstrates that the colour index on the vertical 
axis, (to280 — B), provides a larger spread and more sensitiv- 
ity to the mean age of the stellar population by enclosing the 
4000 A-break. Especially in the presence of noise and scatter 
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Fig. 4. Average K-correction for galaxy types: The mean mag- 
nitude K-correction with redshift is shown as a difference be- 
tween observed frame WFI-R-band and restframe SDSS r-band 
for the average galaxies in the four SED types. 



around the mean colour relation of the sample, type limits de- 
fined on the {B — r)-axis lead to considerable smoothing over 
types defined along the sequence of the relation. It will there- 
fore prove difficult to compare directly our results split by type 
with either 2dFGRS or SDSS data. 

Fig. ^ also highlights a problem with a mono tonic interpre- 
tation of the SED axis among the starburst templates. It seems 
that the templates SB4 and SB5 reside at roughly the same lo- 
cation along the main axis of the distribution, which perhaps 
means physically that their stellar population could have a sim- 
ilar mean age. They differ indeed by their location across the 
main axis, maybe suggesting a second-order feature in the dis- 
tribution of stellar ages. As a matter of speculation, the brighter 
B-band flux combined with a fainter m28o flux might corre- 
spond to a post-starburst galaxy with less current star forma- 
tion as in SB5, but a stronger A star population boosting the 
B-band. However, for the purpose of this paper, we will not fur- 
ther speculate about the detailed star formation histories of the 
templates. We just note, that the starburst range of the templates 
do not form a perfectly monotonic, one-dimensional sequence! 

Furthermore, part of the present sample of galaxies is 
shown split by redshift layer in Fig. ^ depicting the continu- 
ous SED type parameter over luminosity. In this diagram, the 
intervals corresponding to typel to typeA used in the following 
analysis are indicated. A dense horizontal feature can be seen 
around SED values of 85 (=SB4), which is probably caused by 
the not quite monotonic distribution of the starburst templates 
discussed above, which makes a proper statistical SED parame- 
ter estimate very difficult, given that it is based on a probability 
distribution over a locally non-monotonic parameter. 



3.3. Restframe luminosities and sample definition 

Instead of using generic K-corrections, the restframe luminos- 
ity of all galaxies are individually measured from their 17-filter 
spectrum. For each galaxy, three restframe passbands are con- 
sidered, (i) the SDSS r-band, (ii) the Johnson B-band and (iii) 
a synthetic UV continuum band centered at Arcst — 280 nm 
with 40 nm FWHM and a rectangular transmission function. 
This is achieved by precisely matching the redshifted template 
corresponding to the galaxy SED classification into the ob- 
served multi-colour photometry and integrate its spectrum over 
redshifted versions of the three restframe passbands to derive 
the flux to be observed in them. At 0.2 < z < 1.2, this ap- 
proach allows a reliable measurement of the luminosity in the 
B- and 280-band, but does require an extrapolation for the r- 
band at z <; 0.5, where it is redshifted beyond our longest 
wavelength filter The mean K-correction of the average galax- 
ies in each SED type is shown in Fig. ^ as a function of redshift. 
Throughout the paper, we use Hq = h x 100 km/(s Mpc) in 
combination with (fim, = (0.3, 0.7). 

The sample used for all analyses in this paper is defined 
by limits in aperture magnitude, in redshift and in luminosity. 
It is still affected by incompleteness within these limits as dis- 
cussed in the following section. Objects are selected to have 
an aperture magnitude of 17 < R < 24, because in the sat- 
uration range of the individual frames as well as in the noisy 
magnitude regime we can not reliably measure spectral shapes 
and redshifts. They are further selected to have a redshift of 
0.2 < z < 1.2: At low redshift the solid angle covered by 
our survey is too small to obtain useful samples; and at higher 
redshift we currently have no spectroscopic information on our 
redshift accuracy, not even from the CADIS observations of 
^ 100 galaxies. The resulting catalogue contains more than 
25,000 galaxies of which ^ 50 appear to have luminosities of 
Afi < —24, corresponding to 0.2% of the sample. We assume 
that these objects are mostly unreliable measurements, due to 
(i) Seyfert 1 galaxies contaminating the sample with potentially 
entirely wrong redshift estimations, (ii) catastrophic mistakes 
in the redshift assignments leading to completely wrong lumi- 
nosities, (iii) any photometric artifacts of unknown origin. We 
therefore restrict the sample to M,- > —24. For a discussion of 
the average reliability of redshift measurements we like to refer 
the reader to Sect. 3.5. 

3.4. Completeness correction 

The subsequent analysis will draw on galaxy catalogues, con- 
taining only objects with successful z estimates and SED clas- 
sifications. It is therefore critical to understand for which galax- 
ies the data permit such a classification in which fraction of 
cases. 

The subject of this so-called completeness correction of 
a catalogue, or its resulting luminosity function, is indeed 
a fairly complex one: In order to correct a measurement of 
0(-^^tot, z, SED), we would like to have ideally a complete- 
ness function C(Aftot, z, SED) that obviously depends on the 
total absolute magnitude Mtot- 
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Fig. 5. The galaxy sample (here R < 23) split by redshift interval: Shown is the spectral type parameter SED over the luminosity 
Mr in the restframe SDSS r-band for (fim, SIa) = (0.3, 0.7). The intervals occupied by the four spectral types used in this paper 
are delineated by dashed lines. At all redshifts a clear trend between luminosity and SED is demonstrated showing the well- 
known fact that more luminous galaxies are redder on average. The dense horizontal concentration around the SED 85 originates 
from the fact, that the template library is not entirely monotonic (see relative position of templates SB4 and SB5 in Fig. 



On the contrary, the observed survey is characterized by the 
signal-to-noise ratio of the photometry, which determines the 
classification performance and the completeness of the redshift 
catalogue. These all depend on the aperture photometry we use 
to establish the spectral shape, which essentially measures the 
central surface brightness of objects after convolving their ap- 
pearance outside the atmosphere with an effective PSF of 1'.'5 
diameter. 

Given the survey parameters, the completeness of the 
classification and redshift estimation can be derived from 
Monte-Carlo simulations as applied extensively and explained 
in Wolf, Meisenheimer & Roser (2001) and as akeady 
used for the derivation of galaxy luminosity functions in 



CADIS (Fried et al. 2001). The product of these simula- 



tions is a completeness map in fine bins of observed, con- 
volved central surface brightness, redshift and SED type 
(^(/iccntiai, z, SED). If we had proper knowledge of the dis- 
tribution function p(Aftot|Aiccntiai, z, SED), we could derive 



the function C{Mtot, z, SED) needed by a convolution: 

C(Mtot, 2, SED) = C(Mcc„tral, z, SED) * (1) 
p{Mtot\^J■ccntrah S ED) . 

This approach is necessary, because galaxies of lower surface 
brightness which are near but below the total apparent aper- 
ture magnitude limit of the sample will drop out of the sam- 
ple, although formally their overall total magnitude might sug- 
gest inclusion into the sample. Therefore, any completeness 
correction for the luminosity functions will ultimately be con- 
strained by our limited knowledge on the distribution function 
p(Aftot|Mccntiab z, SED), which we try to estimate from the 
sample itself, while being affected by its very incompleteness 
(see Fig. |l§. 

For unresolved galaxies, our classification algorithm is 
more than 90% complete across the redshift range considered 
here for most galaxies with R ^ 24. After convolving the com- 
pleteness map for point sources with the redshift-dependent 
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Fig. 6. Completeness maps for unresolved galaxies: Grey-scale and contour maps demonstrating how the fraction of galaxies 
having successful redshift measurements depends on magnitude, redshift and spectral type. The greyscale shows completeness 
levels from 0% (light grey) to 120% (black). The contour lines are drawn for 90% (white) and 50% completeness (black). This 
completeness map corresponds to C(/iccntrai, z, SED) as in Eq. (1), but plotted over the R-band magnitude of unresolved objects. 
Values above 100% occur when redshift aliasing creates local overdensities in the observed, estimated z/SED-distribution based 
on a flat underlying simulated distribution. These maps are based on Monte-Carlo simulations of the survey. See Sect. 3.4 for a 
detailed discussion of the completeness correction. 



distribution of typical corrections to total galaxy magnitude, 
we find that the sample should still be 90% complete at roughly 
R ^ 23, although this number depends in detail on redshift and 
galaxy type (see Fig. The 50% completeness line ranges 
from R < 23.8 to i? < 23.5 for typel to typeA, respectively. 



3.5. Influence of limited redshift accuracy 

For the subsequent analysis it is crucial to check to which 
extent the limited redshift accuracy (cr^ « 0.03) affects in- 
ferences about luminosity functions, compared to samples of 
high-resolution redshifts obtained from slit spectroscopy. Two 



major aspects need to be explored, redshift aliasing, and catas- 
trophic mistakes, which are both irrelevant for interpreting 
well-exposed data from a spectrograph, but could play a role 
in our case: 

1 . Aliasing results from the presence of structures finer than 
the resolution limit in violation of the sampling theorem, 
and appears as fake structure on a scale that is typically 
slightly larger than the resolution. In this paper, we avoid 
dealing with the problem altogether by not trying to inter- 
pret structure close to the resolution limit (cr^ ^ 0.01 for 
bright typel galaxies and (Tz ^ 0.1 for faint typeA galax- 
ies), and instead choose our redshift bins wide enough. 
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Fig. 7. Completeness maps for resolved galaxies: These full completeness maps for total magnitudes are generated from 
(^(/iccntrai, z, SED) by Selecting only objects of i? < 24 (survey criterion) and convolving the map with the distribution of 
surface brightnesses as found in the survey, resulting in C(i?tot, z, SED). Contour lines are drawn for 90% (white) and 50% 
completeness (black). Obviously, the selection function is fairly uniform within the redshift interval from 0.2 to 1.2 and remains 
on average above 90% for R < 23. The limiting magnitude for 50% completeness is roughly at i? « 23.8 and 23.5 for typel and 
typeA- galaxies, respectively. 



Catastrophic mistakes occur in certain regions of colour 
space where different interpretations can be assigned to the 
same colour vector and probabilistic assumptions are used 
to make a final redshift assignment. Obviously, a number 
of cases will involve an assignment of the wrong redshift, 
but their effect on luminosity functions should be relatively 
small. They would only make a real difference at the most 
luminous end, if many low-redshift objects of medium lu- 
minosity were wrongly assumed to reside at high redshift 
suggesting a high luminosity and consequently boosting the 
abundance of rare giant galaxies. From a simulation we in- 
ferred that the effects will be negligible even if random red- 
shifts were assigned to 10% of all galaxies in the sample. 



At the moment a systematic cross check with significant 
samples of spectroscopic redshifts is still pending for COMBO- 
17. Therefore, we estimate the redshift accuracy drawing on 
the spectroscopic cross-check for CADIS using ~100 galax- 
ies at z < 1.2 (Wolf et al. 2001a), which imply an average 



redshift error of (Tz ~ 0.03. Given a certain similarity of the 
CADIS filter set and similar observing and reduction proce- 
dures, it is natural to assume that the COMBO-17 redshifts 
would be at least as reliable and accurate as the CADIS red- 
shifts. Monte-Carlo simulations of the two filter sets argue in 
fav our of a redshift accuracy that is slightly better in COMB O- 
17 ( |Wolf, Meisenheimer & Roser 200l| ; |Wolf et al. 20"oIa| ). In 
CADIS, we have seen a fraction of catastrophic mistakes at 
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Fig. 8. The full galaxy sample split by field: Shown is redshift over total luminosity in the restframe SDSS r-band for , Ha) = 
(0.3, 0.7). All three fields combined contain25,431 galaxies in the range 0.2 < z < 1.2 (upper left panel). The magnitude limit of 
the sample (R < 24) appears as a fuzzy faint limit on the right side of the distribution. The fuzziness originates from applying the 
magnitude limit not to the total magnitude but to the aperture magnitude of unresolved sources (see also text). Several horizontal 
features can be identified which represent local overdensities. The most conspicuous feature is the concentration just below 
redshift 0.2 in the A901 field, which represents the Abell clusters 901 and 902. Roughly 1000 cluster members are identified in 
the A901 field. The SI 1 field contains another rich cluster at z — 0.11. 



i? < 24 of 10%. However, the classification and redshifts 
in the COMBO- 17 fields are currently checked against spectra 
from 2dFGRS and 2QZ at the bright end, as well as against 
spectra of X-ray sources obtained in the CDFS at the faint end. 
Early indications are, that among bright galaxies at z < 0.3 the 
accuracy is better than Uz < 0.02 and the rate of catastrophic 
mistakes well below 10%. The Seyfert galaxies seen at fainter 
magnitudes in the CDFS are the most challenging objects for 
the redshift estimation, but their rate of redshift outliers is still 
below 20% and their accuracy is still on the order of (Tz ^ 0.05. 

Once, we have obtained better knowledge on the redshift 
error distribution, we should also assess their impact on lumi- 
nosities more carefully. A mean error of 0.03 will lead to less 
than 0™1 scatter in the luminosities at z ~, but up to 0'."25 
scatter at z 0.3. The situation is more dramatic for faint star- 
burst galaxies, where we estimate redshift errors around 0.1, 



implying errors in luminosity of ±0'.°25 at z '--^ 1, increasing to 
±0'."75 at z ~ 0.3. The latter errors potentially bias the steep 
luminosity function of starburst galaxies to brighter L* values. 
However, we reserve a more rigorous treatment of this issue to 
a later publication, following spectroscopic work for establish- 
ing the error distribution. 



4. Results: redshifts and SEDs of ~25,000 
galaxies at 0.2 < z < 1.2 

The immediate result of the above classification and redshift 
determination is shown in Fig. |]and^, where the results from 
the three fields (see Table ^ are combined, and in Fig. ||, where 
the sample is split by field. We restrict our detailed analysis of 
the luminosity function as presented in Sect. 5 to the redshift 
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Fig. 9. The full galaxy sample split by spectral type: Shown is the redshift over the luminosity in the restframe SDSS r- 
band. Several horizontal features represent local overdensities. The most conspicouos feature is the concentration just below 
redshift 0.2, which represents the Abell clusters 901 and 902. The ^^1000 cluster members follow the known relationship of 
lower luminosity for bluer galaxy types. 



interval z=[0.2,1.2], which contains 25,431 measured redshifts 
for galaxies at i? < 24. 

4. 1. SED distribution 

Fig. H shows the sample in terms of SED parameter over rest- 
frame SDSS r-band luminosity, split in four different redshift 
intervals of width Az — 0.2. In this diagram, the magnitude 
limit of the sample appears as a faint limit on the right side of 
the distribution. If a sample was selected on the basis of a total 
SDSS r-band magnitude and observed at z = 0, the selection 
limit would show up as a perfectly vertical border line. 

However, by selecting in the WFI R-band and - more im- 
portantly - by observing redshifted galaxies, the SDSS r-band 
luminosity corresponding to the selection limit depends on the 
restframe colours of the galaxy. Lower values of the SED pa- 
rameter resemble redder galaxies, which are selected in their 
fainter restframe UV at higher redshift and therefore need to be 
more luminous to enter the sample. This explains the angle of 



the faint selection limit to the right side of the galaxy distribu- 
tion. 

Another factor softening the border line is the fact that 
the selection limit applies not to the total magnitude but to 
the aperture flux and ultimately depends on the central surface 
brightness of the galaxy. The selection by aperture flux corre- 
sponds only for unresolved objects to a fixed magnitude value, 
while extended galaxies will have brighter total magnitudes to 
a varying degree that depends on their morphology. The me- 
dian correction for this effect is just below half a magnitude and 
gives the selection limit a fuzzy appearance. Quoting a limit for 
99% completeness strictly speaking requires knowledge about 
the true abundance of low-surface-brightness galaxies, which 
could always have escaped the object detection to some degree 
while featuring high total luminosities distributed over a large 
area. 

The limit on the left side is produced by the steepness of the 
bright end in the luminosity function. Here the figure clearly 
demonstrates the well-known fact, that the highest luminosi- 
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Fig. 10. Distribution of aperture correction: The distribution of 
aperture correction values as a function of redshift. No correc- 
tion applies to unresolved galaxies. A more detailed investi- 
gation addressing issues of central surface brightness, concen- 
tration and size parameters of galaxies will be based on deep 
HST/ACS images covering the entire CDFS field (cycle 12, PI 
Rix) and published in forthcoming papers devoted to the sub- 
ject. 



ties are found among the reddest galaxies, and that there is a 
clear monotonic relationship of L* with SED type. This appar- 
ently smooth relationship suggests, that parametric fits to the 
luminosity functions should be calculated in a bivariate fash- 
ion, depending both on luminosity and SED type or a suitable 
restframe colour. 

4.2. Luminosities and redshifts 

Fig. H shows the sample in terms of redshift over restframe 
SDSS r-band luminosity, split by observed target field. In this 
diagram, the magnitude limit of the sample {R < 24) appears 
again as a faint limit on the right side of the galaxy distribu- 
tion. However, since the total luminosity corresponding to the 
selection limit depends on redshift, restframe colour and mor- 
phology, this border line is again not defined as sharply as an 
aperture flux limit in the observed frame. The different pan- 
els show the imprint of large-scale structure in the individual 
fields. Narrow horizontal stripes point to clusters and sheets. 
While the 3-D positions (x, y, z) of the galaxies can, of course, 
be used to find new clusters out to redshifts of z < 1.2, some 
clusters in the fields have akeady been known. In fact, the A901 
field has been selected for the very reason of studying the clus- 
ters A901/902 in detail, which show up with their '^1000 iden- 
tified cluster members just below z = 0.2. The pencil beams 
of this survey contain a natural mix of environments at all red- 
shifts. In this paper, we do not attempt any differentiation be- 
tween field and cluster galaxies. Clusters do show up at small 



numbers that change with redshift and fluctuate significantly. 
Therefore, our results for typel galaxies could be affected by 
peculiarities introduced by clusters in certain redshift bins. 

Fig. ^ finally shows the sample in terms of redshift over 
restframe SDSS r-band luminosity, split by spectral type. Here, 
the dependence of restframe colour on the SED type affects 
the faint selection limit. At higher redshifts the redder galaxies 
are observed in their fainter restframe UV region and there- 
fore need to be intrinsically more luminous to be included in 
the sample. This explains the flatter angle of the magnitude 
selection for typel as compared to the steeper border line for 
the blue starburst galaxies in typeA. We are able to trace star- 
burst galaxies down to « — 20 all the way out to 2 « 1.2, 
but in contrast galaxies with restframe colours of present-day 
spheroids around z '--^ 1 are only identified at w —21.5. 
The bright cutoff on the left side of the galaxy distribution al- 
ready demonstrates how the evolution with redshift depends on 
the spectral type. Luminous typel galaxies show basically no 
trend with redshift while typeli and 4 galaxies show a strong de- 
pletion of luminous objects when going from z = 1 to z = 0, 
either due to dimming or dropping density. 



4.3. Ttie redsfiift distribution of R< 24 galaxies 

Our subsequent discussion of the galaxies is phrased in terms of 
absolute magnitudes, i.e. luminosities, and redshifts. However, 
for many applications the redshift probability distribution of 
galaxies brighter than some apparent limit is of considerable 
interest. Foremost among such applications is perhaps grav- 
itational lensing, in particular weak lensing, where the small 
image distortions of thousands of galaxies are combined to ex- 
tract information about the intervening mass distribution (e.g. 
Mellier, 1999). In almost all practical cases there is no di- 
rect redshift information available about these numerous faint 
source galaxies. 

In Fig. [n] (left panel) we show the redshift histogram of 
all galaxies from the catalogue discussed above for magnitude 
limits of i? < [22, 23, 24]. These histograms can be corrected 
for incompleteness, using the estimates from Section 3.4 (see 
right panel). The median redshift is z 0.55, with the 90% 
lower and upper percentiles being z — 0.11 and z ~ 1.20, 
repectively. 

Redshift distributions to such faint apparent magnitudes al- 
read exist over very small fields (e.g. the HDF North and South; 
Cohen et al., 2000). However, there the uncertainty in the re- 
sulting histogram is completely dominated by the field-to-field 
variations, whose fractional contrast is near unity in Az ~ 0.2 
bins over such small fields. Here we can present for the first 
time the redshift distribution of galaxies to i? ^ 24 over fields 
large enough that the variance is low. 



4.4. Calculation of the luminosity function 

We follow the usual definition of the luminosity function as 
the number density of galaxies in an interval of luminosity or 
absolute magnitude, expressed in units of h^Mpc^^. Here, we 
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Fig. 11. Cumulative redshift histograms of our galaxy sample: Left panel: Uncorrected cumulative histograms of galaxies found 
in the catalogue at i? < [22,23,24]. Right panel: Redshift histograms corrected for completeness. The median redshift for 
i? < 24 is ^ 0.55 in either sample. In all cases, total magnitudes are used. The correction is based on the completeness maps in 
Fig. 0. Our redshift accuracy at z > 1.2 is currently unknown. 



use two estimators for its calculation (see Willmer 1997 for a 
comprehensive overview): 



4.4.1 . Non-parametric estimates of (/>(M) 

We use the non-parametric 1/Vmax estimator ( [Schmidt 1968 ) 



in the form proposed by Davis & Huchra (1982) and modified 
by a completeness correction as outlined in Fried et al. (2001) 
to calculate the differential luminosity function, which is just 
given by the sum of the density contributions of each individual 
galaxy in the considered luminosity/redshift/SED-type bin: 



(j}{M)dM = 



1 



V^{M,z,SED) 



(2) 



Vi {M, z, SED) is the total comoving volume in which the 
galaxy i could be located to be included in the sample. The 
boundaries of this volume are given by the redshift interval in 
question and by the selection function C{M, z, SED) of the 
survey as discussed in Sect. 3.4. Altogether, the volume is 



V^{M,z,SED) = An 



CiM,z,SED)4^dz 
dildz 



(3) 



The errors of (j){M)dM are just given by the square roots 
of the variances ct^ = yJ^~YJVf\M^7^SWD), which means 
that we take only statistical Poisson noise from the galaxy 
counts into account, but ignore any errors in magnitude or 
redshift here. Furthermore, we ignore the effects of clustering 
here, since the sample is fairly large and should not be strongly 
affected. As we will see later, the Poissonian error bars ap- 
pear to be a fair representation of the scatter in the data points 
around parametric fits of the luminosity function, so we believe 
that we do not underestimate our errors. 

However, due to biases and selection limits the binning pro- 
cess can lead to a suboptimal representation of the underly- 
ing measured luminosity function. As Page & Carrera (2000) 
have pointed out, the faintest bin containing the selection cutoff 



within its bin limits can be significantly harmed. Therefore, we 
decided to ignore any Vniax data points, where the magnitude 
cutoff of the survey shrank the accessible volume of the bin by 
more than 30% compared to an infinitely deep survey. 

4.4.2. Parametric estimates of M* and a 

For a parametric maximum-likelihood fit to a Schechter func- 
tion, we use the STY estimator (Sandage, Tammann & Yahil 
1979). The Schechter function (1976) is 

(j){L)dL = (l)*{L/L*Ye- dL (4) 

or in terms of magnitudes using x = iq^^-M^-m ) 



(j){M)dM = 0.4 In 10 0*x"+ie-"' dM 



(5) 



In the STY formalism, the free parameters of an underlying 
parent distribution (here M* and a) are tested by calculating 
the probability that the observed magnitude distribution is con- 
sistent with the parent distribution (t){M), given the selection 
function C{M, z, SED) of the survey. The probability that a 
single galaxy of luminosity M will be detected is given by 



p{M, z, SED) 



(l){M)C{M, z, SED) 
J (t){M)C{M, z, SED) dM ' 



(6) 



where the integral is effectively limited by the selection 
function containing also upper and lower magnitude cutoffs. 
The joint probability that all galaxies in the sample belong to 
the same parent distribution is then 



C = YIp{M,z,SED) 



(7) 



and this likelihood needs to be maximized by varying M* 
and a. As shown later, we observe an upturn of the luminosity 
function for typel galaxies, just as the 2dFGRS has seen it be- 
fore. We treat this upturn as an extra component and restrict the 
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luminosity range for the STY fit (only for typel galaxies) to ob- 
tain the parameters of the dominant, luminous component, us- 
ing magnitude limits of (M28O, Mb, M^) < (-15, -17, -18). 
At low redshift, the parameters M* and a are well constrained 
and their errors are estimated from the error ellipsoid defined 
as 



ln£ = ln/:„,ax-^X^(^) 



(8) 



where Xp (N) is the P point of the distribution with N 
degrees of freedom (Ax^ = 2.30 and 6.17 for Icr-and 2(T- 
limits corresponding to confidence intervals of 68.3% and 95%, 
repectively). Individual errors on M* and a are quoted by ig- 
noring the covariance but measuring their confidence interval 
within projected contours of Ax^ = 1.0. 

For most restframe passbands and galaxy types we do not 
constrain the knee of the luminosity function out to the highest 
redshifts. The covariance between a and M* makes it virtu- 
ally impossible to obtain a well constrained measurement. We 
therefore decided to measure aiocai at low redshift and assume 
it does not vary with redshift. We then look only at sections of 
the likelihood map in the interval aiocai(type)±(5Q:(type), pro- 
viding us now with a more constrained estimate for M* which 
is valid only under the assumption made, that a does not vary 
with redshift. The error on this M* estimate is obtained after 
rescaling the x^-map such that the best fit M* value found at 
aiocai(type) has Ax^ = 1.0 with respect to the discarded global 
minimum of the full map over the unlimited range in a. The 
error we quote is given by the confidence interval within the 
projected contours of Ax^ = 1.0 above x^(Af*, aiocai) found 
in the map constrained in a. 

4.4.3. Determination of (p* 

Since the normalisation 0* of the luminosity function cancels 
in our implementation of the STY fit procedure and only the 
shape of the distribution is subject to the test, we calculate (p* 
afterwards by 



(9) 



/^f 0(M)/0* dM 



where the average galaxy density n within a certain volume 
is given by 



»gal 



(10) 



There are two different error sources for 0*, one arising 
from cosmic variance, i.e. large-scale structure, and the other 
from the covariance with M* . All errors quoted for (/)* in this 
paper are determined from field-to-field variation among our 
three fields, whatever redshift bin or SED type we consider; this 
turns out to be the dominant source of uncertainty. However, 
we do list the covariance with M* as well, but note, that the 
latter does not dominate over the error from structure under the 
assumption made that a does not vary with redshift. 



4.4.4. Determination of j 

We furthermore calculate the luminosity density of galaxies, j, 
after integration over the luminosity axis. The full integrals can 
easily be obtained from the STY fits as, e.g.. 



j = y dL<P{L)L = 0*L*r(a + 2) 



(11) 



where V is the Gamma-function. Later we show explicit 
results for a range of finite lower luminosity limits to illus- 
trate the effect of extrapolating the luminosity function with 
constant a for the different types. The very assumption of a 
remaining constant to infinitely low lumniosities is the largest 
source of error for j -estimates at high redshift. It is not reflected 
in any error estimates we quote, but can be assessed from the 
convergence of the integral with decreasing luminosity limits 
as shown in the following section. If we believe the STY fits 
and use their formal errors to derive an error estimate on j, we 
find the error in cf)* originating from field-to-field variation to 
dominate and derive j errors directly from these. As solar lumi- 
nosities we use (Vega-normalised) values of Mjgg = 6.66, 



B,0 



5.30 and NL 







4.47. 



5. Discussion 

5. 1. The "quasi-local" luminosity function atzr^ 0.3 

With a total field size of 0.78 our survey volume at z < 0.2 
is too small for a sensible derivation of a luminosity function. 
Also, the present redshift errors can have a larger impact onto 
the luminosity measurement in this very local regime. Our low- 
est redshift bin therefore covers the range z = [0.2,0.4] and 
contains 5,674 galaxies from three target fields. Already the 
median redshift in this bin is 0.34 in the luminosity range of 
A/r — [—23, —17]. This value is higher than for recent, large 
surveys of the "local", or present-day, galaxy population, such 
as SDSS or 2dFGRS. We therefore refer to our lowest z, refer- 
ence bin as the quasi-local sample. 



5.1 .1 . Tlie field-to-field variation of the luminosity 
function 

In our lowest redshift bin, z = [0.2,0.4], the co-moving vol- 
umes are smallest and the structure in the galaxy popula- 
tion has developed furthest. Therefore, this redshift bin is the 
most critical to explore to which extent our inferences about 
the galaxy distribution are affected by field-to-field variations. 
Fig. shows the luminosity functions in the restframe SDSS 
r-passband derived separately from the three individual fields 
and for the whole sample combined. All error bars shown re- 
flect only the la Poisson variance from the finite number of 
galaxies in the bin. It is apparent that due to large scale struc- 
ture the results differ somewhat more from field to field than 
implied by Poisson errors. However, for the most part, these 
differences are merely fluctuations in (j)*, which we use to esti- 
mate errors on ip* . 

The CDFS field appears to have a lower galaxy density than 
the other two fields. We note, that observations of the Chandra 
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Fig. 12. Field-to-field variation: Comparison of the luminosity function (/'(Mi ) of galaxies from the lowest COMBO-17 redshift 
bin (z = [0.2, 0.4]) among the three disjoint survey fields. Error bars are 1-cr-Poissonian. Grey line: STY fit of Schechter function 
to the whole sample. Black line: fit to the panel-specific field sample. By choice, some galaxies in the A901 field should be subject 
to magnification by gravitational lensing. We estimate, that this brightening corresponds to Am ~ — 0™25 in the central density 
peak of the cluster, but less in its outer regions. At this point, we neglect the effect for the field as a whole and assume that is 
does not alter the conclusions of this paper The hump around il/ ^ — 21 in the LF of the A901 field appears to originate from 
large-scale structure, as it is related to a group of bright galaxies at similar redshift near the Eastern edge of the field. 



observatory have revealed that the CDFS also shows number 
counts of X-ray sources lower b y a factor of two compared to 
the CDF-North ( [Norman 2002| ). The implication of this field- 
to-field variance is complicated as the fields were not chosen 
at random: while the CDFS is very "empty", the A901 field 
was chosen to contain a cluster (though not in the redshift bin 
shown). This test shows, that even with a 0.5° x 0.5° field size 
the field-to-field variations are still noticeable and must be ac- 
counted for, but they are much less pronounced than in tiny 
fields, such as the Hubble Deep Fields. 



5.1 .2. The SED type and passband dependence of 
the luminosity function 

The ^5,700 galaxies in the quasi-local bin are enough to draw 
up a quite comprehensive statistical picture of the (z) « 0.3 
galaxy population, which is summarized in Fig. |l3[ We have 
enough galaxies that we can study the LF in the four broad 
SED classes defined in Section 3.2. We have enough wave- 
length coverage that we can construct the LF in three wave- 
bands, 280 nm, hj, and r. Further, with i? <^ 24 our data reach 
well below the knee in the LF, M*, so that the faint end slope 
a is well constrained. 

Is is obvious from Fig. |l3| how much the faint end of the 
luminosity function depends on SED type: the later or bluer the 
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Fig. 13. SED dependence of the luminosity function: Results are shown for quasi-local sample at z = [0.2, 0.4] in three different 
restframe passbands, SDSS-r, B, and at 280 nm, from top to bottom respectively, and split by SED types. Left panels: Vmax 
datapoints with Poissonian error bars and STY fits of Schechter functions. Right panels: Likelihood contours for parameters M* 
and a with la and 2a error contours. Differences in M* between the bands reflect mean restframe colours of the type samples. 
The upturning faint end of typel has been excluded from the STY fits. Note the relatively fainter M* of typeA (starburst) galaxies 
compared to type3. 
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Fig. 14. COMBO-17 vs. SDSS: the luminosity function (t){M^) of galaxies from the quasi-local COMBO-17 sample at z 0.3 
compared with the SDSS results by Blanton et al. (2001), which are represented as error bars connected by a dashed line. The 
original AB magnitudes are converted to Vega magnitudes, and the types have been matched by adjusting the SDSS {g — r )-colour 
cuts, although this will not lead to an emulation of our type definition (see text). For the combined sample of all SED types, the 
basic match around the knee of the function is quite good, but differences appear at the bright and faint end. While faint-end 
differences relate to evolution of faint blue galaxies between the different median redshifts of the two faint end samples, ^ 0.34 
and ~ 0.05, bright-end differences might result from different sensitivity to faint surface brightness envelopes in giant spheroids, 
as they are expected between 5 -hour and 1 -minute exposures. 



SED type, the steeper the faint end of the luminosity function. 
This effect, quantified in the right hand panels of Fig. |l3| is 
present in all three wavebands. We remind the reader that the 
SED types are defined by redshift- and luminosity-independent 
restframe colour, not by morphological or evolutionary type. 

For the most part, these SED types reflect a mean stellar age 
sequence, although metallicity effects complicate the relation- 
ship between age and colour: in particular for early types the 
colour-luminosity relation will place fainter (presumably more 
metal-poor and hence bluer) galaxies of a given stellar age pref- 
erentially into a later SED bin. This effect can contribute to a 
downturn at the faint end of type 1 , but it will by far not be the 
dominant effect leading to the positive a. 

Not surprisingly, the characteristic absolute magnitude, 
M*, of the LF depends both on SED type and on the ob- 
served waveband. In the r-band M* is nearly independent of 
the SED for typel to 3, while the starburst galaxies of typeA 
are significantly fainter. When looking at the other bands, the 
mean restframe colours of the galaxies belonging to the dif- 
ferent types (as shown in Fig. ^ naturally shift the M* values 
such that bluer galaxies appear relatively more luminous in the 
bluer bands. In the hj band and at 280 nm, the brightest M* 
is found for type3 objects and probably originating form large, 
strongly star-forming spiral galaxies. 



Obviously, we do not observe the entire galaxy population, 
but only galaxies with M Ss —16, i.e. galaxies to SMC lu- 
minosity, in the quasi-local sample. While our parametric fits 
and forthcoming density calculations assume the LF to con- 
tinue to infinitely faint levels at constant a, the LF might turn 
over at some lower luminosity or simply have a more com- 
plicated shape as given by the Schechter function, especially 
when we look at the combined sample of all SED types. The 
sum of Schechter functions for individual types with different 
(M*, (/)*, a) -parameters will not in general give a Schechter 
function again, unless the steepest function dominates the sum 
everywhere as it is almost the case in the 280-band. In Fig. [T^ 
we can see the bright end of the B-band and r-band LFs be- 
ing dominated by typel and the faint end by typeA. At inter- 
mediate luminosities we get almost equal contributions from 
all types and the result does not resemble a Schechter function 
so clearly anymore. Indeed we can see a variation of a with 
the luminosity domain since the typeA LF contributes signifi- 
cantly only in the faintest domain. If typeA galaxies continue 
to rise at their rate fainter than Af280 ^ —16, they will drive 
up the a-value for the combined sample as well. The very rea- 
son for our observed a-values, which are steeper than those for 
the 2dFGRS- and SDSS-samples (see following section), is the 
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Fig. 15. COMBO-17 vs. 2dFGRS: the luminosity function 0(Mb) of galaxies from the quasi-local COMBO-17 sample compared 
with the 2dFGRS results by Madgwick et al. (2002), which are transformed into the Johnson-B system and represented by a 
dashed line (error bars are vanishing in size). The original 2dF galaxy types have been adopted without change and therefore 
correspond only very roughly to our types. The agreement of the combined sample is rather good with slight differences at the 
faint end arising from evolution of faint-blue galaxies. 



stronger prominence of starburst galaxies at the higher redshifts 
we look at. 



5.1 .3. Comparison of COMBO-1 7 with 2dFGRS 
and SDSS 

Figures |l^ and |l5| present a comparison of the COMBO-17 
quasi-local sample to the yet more nearby samples from SDSS 
dBlanton et al. 200 1| ) and 2dFGRS ( [Vladgwick et al. 2002[ ), 
which both have {z) « 0.1. For the 2dFGRS we plot the lu- 
minosity function just for their published types without any at- 
tempt to adjust the respective SED type definitions. Therefore, 
differences are expected to some degree. For the SDSS, we 
have obtained the values of the luminosity function in fine bins 
over the [g — r)-colour axis, and have chosen type limits with 
the aim of matching the shapes of their LF to ours. However, 
a brief look at Fig. ^ should remind the reader that cutting the 
sample by intervals on the {B — 7-)-colour axis does not allow 
to tune the limits such that truely similar type samples are ob- 
tained as with our own definition. 



When taking all SED types (leftmost panel), the agree- 
ment with SDSS and 2dFGRS is quite good around the knee 
of the Schechter function. The y„iax results from COMBO- 
17 are shown as squares with errors bars having no ticks, and 
the STY fits as solid lines following a Schechter function. The 
comparison surveys are shown as another solid line connect- 
ing their Vmax data points. For the SDSS, the individual Vmax 



data points are further shown as error bars with ticks, but for 
the 2dFGRS the error bars are omitted because their sample of 
~ 100,000 galaxies renders them virtually invisible. 

The luminosity axis has been adjusted for both surveys, by 
converting from ABmag to Vega-mag in case of the SDSS, and 
by transforming from the photographic fej system to Johnson- 
B in case of the 2dFGRS. Using synthetic photometry on the 
galaxy templates, we derived a mean magnitude offset B — bj 
of 0.12, 0.11, 0.07, 0.06 and 0.09 for typel, type!, type3, typeA 
and the combined sample, respectively. 

There are significant differences to SDSS and 2dFGRS, for 
very faint and for very bright galaxies. These differences could 
be physical or results of the data treatment, and we cannot give 
conclusive explanations. However, two plausible reasons come 
to mind: the mean redshift of dwarf galaxies in the SDSS and 
2dFGRS sample is z « 0.05, while it is z « 0.34 in the quasi- 
local COMBO-17 sample. I.e. there is a difference of over two 
Gigayears in the mean look-back time. Given the strong red- 
shift evolution of dwarf galaxies (see also Sect. 5.2), this dif- 
ference reflects most likely the physical evolution of faint star- 
forming galaxies found by Ellis et al. (1996) and Heyl et al. 
(1997). 

For the brightest galaxies, in large part early type galaxies, 
there may be a completely different explanation for the dif- 
ference: all magnitudes involved are ultimately some version 
of isophotal magnitudes. The luminosities in COMBO-17 are 
based on the multi-hour R-band exposure, which reaches much 
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fainter surface brightness limits (the SDSS exposure is shorter 
than 1 min), and may therefore measure a larger fraction of the 
low-surface brightness, outer parts of early type galaxies. 

While these effects might have caused differences between 
COMBO-17 on one side and 2dFGRS as well as SDSS on 
the other side, we would not conclude to see major differ- 
ences between SDSS and 2dFGRS as their deviation from the 
COMBO-17 LF at z ^ 0.3 is quite similar Early claims of 
significant differences between SDSS and 2dFGRS (Blanton et 
al. 2001) have since been reconciled by Norberg at al. (2002), 
who demonstrated that applying necessary colour transforma- 
tions and further data adjustments with due diligence produces 
compatible results. 

5.2. The redshift evolution of the luminosity function 
since z 1.2 

In this section we extend the discussion to the redshift inter- 
val z=[0.2,1.2], which contains 25,431 measured redshifts for 
galaxies at i? < 24. In this section, we show only STY fits 
(Fig. |l^ for the luminosity functions, and use them as a basis 
for our discussion since they form a good representation of the 
VJnax data points. For reference, we show the full set of plots 
depicting Knax data points together with STY fits for all red- 
shifts, all SED types and all restframe bands into an Appendix, 
where the interested reader can look at finer details and find 
tables with all parameters. 

Luminosity function fits to different redshift bins of the 
same SED type, where we allowed M* and a to vary freely, 
showed that the faint end slope a does not significantly evolve 
with cosmic epoch. We therefore adopted the hypothesis that 
the faint end slope only depends on SED type, but not on cos- 
mic epoch. This approach avoids the ill-defined covariances 
with the M* fit, as less and less of the sub-L* regime is fit to- 
wards higher redshifts. In practice, the STY fits were obtained 
by measuring a in the quasi-local sample and keeping it fixed 
for all other redshifts. While there seems to be an indication of 
a flatter a at z = [0.4, 0.6], the quasi-local a again fits better 
at z > 0.6 (see Appendix). At this point we suspect that the 
apparent, slight flattening of a in the epoch around to z « 0.5 
originates from an unknown and still uncorrected feature of our 
data. 

Using a fixed a for all redshifts works quite well for the 
luminosity functions of individual SED types. But due to the 
change in the type mix with redshift, this approach is useless 
for the combined sample of all types. There, a has been de- 
termined in each redshift interval independently. But at higher 
redshift the flatter low-luminosity regime in the LF is less con- 
strained by the observations and a correct determination of a is 
difficult. 

The evolution of the galaxy luminosity function is best ob- 
served by type. The four different SED types show different 
evolutionary patterns although there are trends linking them. 
Given that the Schechter function contains just three free pa- 
rameters, (M* , (f>* , a), keeping a constant with redshift means, 
that we can observe the evolution of the LF as an (A/*, (j>*)- 
vector evolving with redshift z (see Fig. p7[). 



If any colour evolution lead only to migration of objects be- 
tween types, but left the distribution of restframe colours within 
an SED type sample unchanged with redshift, it would suffice 
to look at only one restframe passband. The other two would 
then contain purely redundant information, because our SED 
types are defined by restframe colour. However, due to finite 
width of the colour intervals associated with the types, slow mi- 
gration can shift the mean colour within types, and still cause 
slightly different results among the three passbands. In prin- 
ciple, such a migration could also lead to a shift in a, which 
again motivates a bivariate luminosity function depending on 
restframe colour directly. 



As shown in Fig. 17, all types have some luminosity fading 



with cosmic time in common but show quite different trends in 
density evolution. For typel, we see a consistent <; 10-fold in- 
crease in density for all three bands when going from z — 1.1 
to z = 0.3. The steep increase suggests that we have not yet 
seen the epoch of maximum density for typel galaxies, and 
will see a continued increase in the future of the universe. 
Furthermore, all three bands show evidence of strong fading 
— from Am « 1.0 mag in the r-band to Am « 1.5 mag at 
280 nm across the redshift range from z ^ 1.1 to z ^ 0.3. 
The stronger fading in the UV is basically a consequence of 
the mean colour in the typel sample getting redder with time 
as the member galaxies age further, a trend also reflected in the 
drift of the mean template for typel objects with redshift (see 
Fig.|). 

For typel, we again find uniform evidence from all three 
bands for fading at the ^ 1 mag level with no indication of 
a change in mean restframe colour The density trend of typel 
galaxies is positive at higher redshift and turns slightly negative 
towards more recent epochs, suggesting an enclosed epoch of 
maximum density for typel galaxies around z « 0.5 . . . 0.7. 

The luminosity function of type3 galaxies shows basically 
no evolution at z > 0.5, but a later reduction in density of 
^ 0.25 dex, meaning that almost half of the galaxies have dis- 
appeared from the type3 bin by z 0.3. 

The starburst galaxies of typeA show a stronger decrease 
in luminosity, ranging from Am w 1.0 mag in the r-band to 
Am « 1.5 mag at 280 nm across the redshift range from 
z ^ 1.1 to z ^ 0.3. The stronger fading in the UV is again 
a consequence of the mean colour in typeA getting redder with 
time. This is not reflected in any change of the mean template 
with redshift, probably due to the non-monotonic behaviour of 
the starburst templates. The trend in colour is consistent with 
the starbursts getting relatively less prominent within the given 
underlying galaxy. Alternatively, one could speculate about in- 
creased dust reddening. The typeA galaxies show no strong 
trend in density in any band. The steep faint end of the typeA 
luminosity function leads to a strong covariance between M* 
and (/)* and an increased difficulty for perfectly disentangling 
density evolution from fading. However, in terms of total lumi- 
nosity, the epoch of maximum starburst activity certainly oc- 
cured somewhere beyond our redshift limit of z = 1.2. 

On the whole, it appears that these patterns can be under- 
stood as a decrease in starburst activity with time following an 
activity maximum at z <; 1 beyond the limits of this work, 
in combination with a propagation of the galaxies through the 
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Fig. 16. Redshift evolution of the luminosity functions ), (t){B) and (j){M2so)- STY fits of Schechter functions in five redshift 
intervals of width Az = 0.2, centered at z = [1.1, 0.9, 0.7, 0.5, 0.3]. The first and last of these fits is plotted in black, and the 
three enclosed intervals are shown in grey. For the individual SED types the faint-end slope a was determined in the "quasi-local" 



sample and kept constant for all redshifts (for a values see Appendix, Tab. 



and A..4). 
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Fig. 17. Redshift evolution of the parameters M* and 0* : Results of STY fits in five redshift intervals centered at z = 
[1.1, 0.9, 0.7, 0.5, 0.3]. For individual SED types the faint-end slope a was determined in the quasi-local sample and fixed for aU 
redshifts. Plotted are M^^^^ — 5 log h values and l-cr error contours derived from the error in M*, the field-to-field variation in 
0* and the covariance between M* and cj)*. The most distant sample (z — 1.1) and the quasi-local sample (z = 0.3) are plotted 
in black, and the redshifts in between in grey. For orientation we note, that in every panel the quasi-local sample is found at the 
faintest (right-most) M* position. 



types as the mean age of the stellar population increases. As 
the starburst activity continues to drop, the fraction of galaxies 
with only old populations continuously increases. 



The comparison with local surveys from the previous sub- 
section suggests a simple continuation of the trends found be- 
tween (z) — 1.1 and (z) — 0.3, such that the density of type I 
galaxies does continue to rise and typeA galaxies continue to 
vanish or fade. Type! and 3 galaxies remain almost unchanged, 
but we remind the reader that our comparison with local sur- 
veys is only of limited value due to significant differences in 
the type definition. 



If the differentiation between fading and density trends is 
too unreliable in places, given only limited constraints on the 
knee of the Schechter function at higher redshift, the way out 
is to look at figures of integrated luminosity density, which we 
do in the following section, where we continue to include the 
local samples into the discussion. 



5.3. Evolution of the co-moving luminosity density 
at different rest-wavelengths since z ~ 1.2 

As an alternative to discussing luminosity functions, we can de- 
scribe the redshift evolution of the co-moving luminosity den- 
sity stemming from galaxies of different SEDs. I.e. we can ex- 
plore the evolution of the integral over the luminosity function, 
which avoids the problems arising from the 0*, M* and a co- 
variances. 

Conceptually, this approach suffers from the fact that the 
faint end of the luminosity function eventually becomes unob- 
servable for a given redshift, and that an extrapolation for this 
unobserved galaxy contribution becomes necessary. However, 
in the present context this is not a serious limitation: For one, 
our data show directly that the faint end slope covers — 1.5 < 
a < 0.5, depending on SED type, so the LF integral will con- 
verge quite rapidly for luminosities below M*. Further, our 
data reach considerably fainter than previous deep surveys, so 
for any given a much less extrapolation is necessary. 

Fig. 18 shows that the integrated luminosity density in the 



different wavebands is quite well constrained by direct observa- 
tions. Five different lower luminosity limits have been applied 
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for the integration, ranging from M < — 18 to Af < —10 to 
show the effect of the extrapolation for different types. In any 
bandpass the luminosity density integral has largely converged 
if the integration interval extends to il/ <^ — 16. 

In all three bands we see roughly the same picture for the 
individual types but significant differences for the combined lu- 
minosity of all galaxies, as the mix of types changes with red- 
shift. As i oc L*(t>* , the former statement is just a consequence 
of the observation, that already the evolution of M* and was 
consistently seen among the different bands. The B-band dis- 
plays most clearly the overall trend in luminosity density from 
z ^ 1 to now: a decrease among strongly star-forming galax- 
ies is accompanied by perhaps a maximum of type! galaxies at 
intermediate redshifts and a ten-fold increase in typel objects. 

The overall evolution of the luminosity density is summa- 
rized in Fig. |9[ showing our own results (large filled symbols) 
in comparison to earlier work (open symbols) and large, con- 
temporary surveys of the local universe (small filled symbols). 
The values of older work have, of course, been adjusted to the 
cosmological parameters used here. 

In the B- and r-band the luminosity density remains virtu- 
ally unchanged from z = 1.1 to z = 0.5, and even j28o drops 
only httle down to z = 0.5. All bands show a subsequent de- 
crease in luminosity density to the present epoch, most strongly 
at 280 nm. Between z ~ and z ~ 1 we observe j280 to drop 
by a factor of ~ 6, while jb and drop only little. If we con- 
sidered our data at z = 0.5 to be outliers towards the top, the 
flat domain would reduce to z = 0.7 ... 1.1 in the B- and r- 
band. The 280-band would then show a more gradual change 
in the slope from the flatter domain around z ~ 1 to the steeper 
gradient in the local universe. 

We can not quite confirm the steep increase oberved by 
Lilly et al. (1996), that leads to apparently inconsistent data 
at high redshift. By using the Loveday (1992) as a local ref- 
erence, their data suggest an increase by a factor of 4.6 from 
z = to z ^ 1 in jB, which is much stronger than what we 
see. However, we believe the results are consistent within their 
error bars. 

Finally, we show what fraction of the luminosity density at 
a given wavelength arises from which SED type. This is illus- 
trated in Fig. |o|for all three wavelengths. While the strongly 
star-forming galaxies of type3+4 produced the bulk of the ra- 
diation in all three bands at z 1, they nowadays contribute 
only an important fraction of the luminosity when considering 
the near-UV (280 nm). On the contrary, the fraction of light 
arising from early-type SEDs {typel) has increased by large 
factors in each of the bands. 

5.4. Comparison with ottier tiigiier-redstiift surveys 

In this section, we would like to compare the findings on 
the type-dependent evolution of galaxies at z < 1.2 with 
other recent surveys, i.e. the Canadian French Redshift Survey, 
CFRS ( pfly et al. 1995D , the CNOC2 s urvey of the Can adian 
Network for Observational Cosmology ( |^in et al. 1999 ), and 
the Calar Alto Deep Imaging Survey, CADIS (Fried et al. 
2001). Our comparison will be done qualitatively, taking into 



Table 4. Characteristics of galaxy redshift surveys compared 
here. Note: We ignore the local CFRS sample at z < 0.2 due to 
its small size, although results have been published. The most 
distant bin of its red sample is actually at median redshift ^ 
0.87. 



Survey 


Median 


of bins 


objects 


Types 


a(z) 




min. . 


.max 








CFRS 


~ 0.35 


. . . 1.1 


~ 600 


2 




CNOC2 


~ 0.2. 


.0.45 


~ 2, 200 


3 


const 


CADIS 


~ 0.4. 


..0.9 


~ 2, 800 


3 




COMBO-17 


~ 0.33 


. . . 1.1 


~ 25,000 


4 


const 



account the differences arising from incompatible type defini- 
tions. 

We summarize beforehand, that altogether the observations 
of all four surveys are consistent with each other and no signif- 
icant contradictions have been seen exceeding the size of the 
error bars. However, COMBO-17 constrains the evolution of 
galaxies much more clearly than any of the comparison sur- 
veys, due to larger sample size and/or wider redshift cover- 
age. Especially, the smaller samples of previous surveys did 
not allow to see any second-order effects as demonstrated by 
the peak in typel density. Instead they had to rely on fitting 
first-order models with redshift to the data. 

Whenever seemingly contrary results were derived, their 
origin is mostly to be found in different type definitions or as- 
sumptions on parameter evolution. We therefore believe, that it 
is not worth the effort to reproduce all the different type defini- 
tions between the surveys in order to do a precise quantitative 
comparison, and omit this exercise. 

5.4.1 . Characteristics of comparison surveys 

All three comparison surveys derived the luminosity function 
in the restframe B-band. Most relevant for the following dis- 
cussion are these characteristics (see Tab. 

- The CFRS covers the same wide redshift range as 
COMBO-17, z = [0.2, 1.2], and features the same median 
redshift of its sample ((z) = 0.56), but provides the small- 
est of the samples considered here with 591 objects. The 
results are derived for two galaxy types with a — a{z), 
a red sample corresponding roughly to our typel+2, and a 
blue sample resembling our type3+4. 

- The CNOC2 reaches only out to z = 0.55, the median 
reshift of COMBO-17 and the CFRS, and also less deep, 
being limited to Rc < 21.5. However, it contains over 
2000 galaxies and has derived results with fairly small error 
bars and a = const for three types of galaxies. Assuming 
that a depends monotonically on restframe colour, and 
comparing the values found for the LF by type, we could 
very crudely identify their three types with our type"l-lS', 
type"l.5-2", and type3+4. 

- Finally, CADIS covers almost the full redshift range of 
COMBO-17 and CFRS, although only three redshift bins 
were used with median redshifts ranging from 0.4 to 0.9 as 
opposed to 0.3 to 1.1 for COMBO-17. Like CNOC2 it con- 
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Fig. 18. Redshift evolution in luminosity density j,-, je and .7280: j integrals calculated from STY fits. The five lines belong to 
different values of the integration limit of M < [—18, —16, —14, —12, —10] and converge to a total luminosity for every value of 
a > —2 (see Eq. 1 1). For typel and jb have already converged at M < —18 rendering the five lines indistinguishable. With 
a « —1.5 for typeA galaxies, the integral to the faintest levels contains up to 2.5 x the luminosity measured to M < —18. Local 
data from SDSS (r-band), 2dFGRS (transformed into B-band) and from Lilly et al. (1996) as based on Loveday et al. (1992) in 
the 280-band. 
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Fig. 19. Overall redshift evolution of the luminosity density in all bands: The data points represent the total luminosity density 
(integrated to faintest luminosities) for 280 nm (circles), B (squares) and r (pentagons). The large solid symbols are the values 
from this paper, but split into finer redshift bins of Az = 0.1 for clarity, the small solid symbols at the low redshifts represent 
the local comparison surveys (Blanton et al. 2001; Madgwick et al. 2002), and small open symbols represent the results from 
previous redshift surveys (Madau, Pozzetti and Dickinson 1998 listed values for Ellis et al. 1996, Connolly et al. 1997, and 
Lilly et al. 1996 also listed values for Loveday et al. 1992). Luminosities of M\ — objects for conversion into physical units: 
L, = 3.84 X 1013 w/Hz, Lb = 4.93 x lO^^ W/Hz andisso = 1-18 x lO^^ W/Hz. 
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Fig. 20. Redshift evolution in luminosity density fraction by type: Fraction of the total luminosity density contributed by a single 
SED type (typel at the bottom, typeA at the top). The total luminosity integrals are calculated from the STY fits. Local data: 
2dFGRS and SDSS. 



tains more than 2000 galaxies, and results were obtained 
with a — a{z) for three types, corresponding roughly to 
our typel, typel and type2>+A. 



5.4.2. Redshift evolution of a? 

The CNOC2 survey considers a to be not changing with red- 
shift, but the CFRS and CADIS allow a to change while inves- 
tigating a much wider redshift range than CNOC2. Although 
COMBO-17 studies a wide redshift range as well, we found 
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that the assumption of a being constant in redshift is reason- 
able, provided that the types are not too broadly defined. 

Indeed, we are compelled to believe, that a is a monotonic 
function of restframe colour, but does not change with redshift 
at least at 2; ^ 1. Any galaxy type encompassing a wide range 
of colours can obviously experience a strong evolution in its 
sub-type mix. In that case, the a of the whole type sample is 
dominated by different sub-types at different epochs and there- 
fore reflects their respective, different a values. 

We like to stress again, that if the steep LF of starburst 
galaxies does not turn over at some low luminosity, it will al- 
ways drive the a of any whole sample to a steep value at the 
faintest end, which should even be observable in local samples 
if only galaxies of sufficiently low luminosity were observed. 

It is this assumption of a being constant, combined with a 
fairly dense sample, that allows the CNOC2 and COMBO-17 
to derive the clearest evolutionary patterns. All three previous 
surveys and possibly even COMBO-17 might have too small 
and shallow samples to measure positively any evolution in a. 
Allowing it to vary freely in CFRS and CADIS has then greatly 
reduced the constraints on the evolutionary patterns they de- 
rived. 

5.4.3. CFRS vs. COMBO-17 

When looking towards high redshift, the CFRS reports to see: 

- no evolution of the luminosity function (neither in ^* nor 
in M*) in its red sample 

- a significant amount of brightening and steepening in the 
luminosity function of the blue sub-sample. 

Their red smaple is at all redshifts dominated by our typel 
galaxies which show only little evolution that could easily be 
invisible in a small sample. The fading with time observed by 
us might have been compensated by an increasing density of 
luminous typeX objects, rendering any changes on the luminous 
side and around the knee of the combined typel+2 luminosity 
function hard to detect. 

Their blue sample is dominated by our type3 galaxies at 
Mb < -20 and by our type4 galaxies at Mb > -20. When 
looking at the steep, luminous end of the LF the density in- 
crease of type3 objects observed by us appears like a brighten- 
ing as well. At the faint end, the brightening of type4 objects 
seen by us causes an increasing fraction (with redshift) of type4 
objects being seen in the faint domain of the combined type3+4 
luminosity function. This is an illustrative case of a change in 
sub-type mix mimicking a steepening. 

5.4.4. CN0C2 vs. CFRS and COMBO-1 7 

The CNOC2 concluded to see (when looking back in redshift): 

- strong rise in luminosity but a reduction in density for its 
early-type sample 

- gentle rise in luminosity but no trend in density for its in- 
termediate sample 

- no trend in luminosity but a strong rise in density for its 
late-type sample 



They fit an evolutionary model to their samples whereby 
M*{z) - M*(0) = -Qz and logp(^)//9(0) = OAPz with 
type-dependent evolutionary parameters P and Q. They derive 
values supporting the clear evolutionary pattern from above, 
albeit with large error bars due to the short redshift baseline. 
In principle, we can predict changes M*(l.l) — M*(0.3) and 
log p(l.l)/ p(0.3) with respective errors, as we should see them 
in COMBO-17. After trying to adjust the mismatch in type defi- 
nitions, the results become rather consistent within large errors. 

For its blue sample CNOC2 found a clear </>* drop with little 
M* change, while the starburst galaxies in COMBO-17 appear 
to evolve only in M*, but there is indeed no contradiction, as 
the CNOC2 blue sample is dominated by our type3 objects, 
which are more prominent at low redshift and < 21.5. 

The apparent differences between the CNOC2 and CFRS 
results are best explained by the difference in type definition. 
For the red CFRS sample, also a partial cancellation of oppos- 
ing trends in density and luminosity plays a role. When com- 
bining part of the intermediate CNOC2 sample with its late- 
type sample, it is easy to imagine that on the whole brightening 
and steepening is seen, again. 

5.4.5. CADIS vs. CFRS, CN0C2 and COMBO-17 

CADIS concluded to see (when looking back in redshift): 

- no trend in luminosity but a reduction in density for its red 
sample 

- no trend in luminosity and only a sUght reduction in density 
for its intermediate sample 

- a significant amount of brightening and steepening in its 
blue sample. 

The CADIS results have been claimed to be mostly consis- 
tent with the CFRS results while improving on accuracy due 
to an increased sample size. They appear to be in disagreement 
with CNOC2, involving a dispute on the justification of keep- 
ing a fixed with redshift. We believe, however, that the obser- 
vations are consistent, and the conclusions could become com- 
parable after matching up the type definitions and agreeing on 
the behaviour of a. Keeping a free also makes it more difficult 
to separate trends in M* and </>* . 

Again, the brightening and steepening observed for the blue 
sample resembles the combined evolution of type3+4 with the 
a value of type4 being more important at higher redshift and 
that of type3 being dominant at low redshift — at least within 
the luminosity interval observed. We repeat here our claim that 
the COMBO-17 observations do not support any change of a 
with redshift for a sample with any given narrow interval in 
restframe colour. It is simply changes in type mix occuring in 
broadly defined types that lead to the observation of a changing 
in the typical luminous domjiins as they are currently observed 
in high-redshift surveys. 

6. Summary and Conclusions 

Based on photometry in 17 (mostly medium-band) filters ob- 
tained by the COMBO- 1 7 survey on three independent fields of 
0.26 each, we have derived redshifts and SED classifications 
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for ~ 25, 000 galaxies to 24. The redshifts may be viewed 
either as high-precision photometric redshift estimates or very 
low resolution (R^ 12) spectral redshifts. They have a preci- 
sion of cr^ « 0.03 and lie mostly in the range 0.2 < z < 1.2 
with a median redshift of 0.55. The upper Umit of galaxy red- 
shifts considered here is set by the availability of spectroscopic 
cross-checks confirming their precision and accuracy. Our re- 
sults are unaffected by random redshift errors of ~ 0.03, and 
are robust even for somewhat larger, systematic errors. 

Compared to previously pubhshed large-sample surveys 
with better redshift accuracy, COMBO-17's flux Umit is nearly 
two magnitudes fainter and the number of galaxies is more 
than an order of magnitude larger. Taken together, the prop- 
erties of the survey allow us to draw up a comprehensive pic- 
ture of how the luminosity function and the SEDs of galaxies 
have evolved over the last half of cosmic evolution. The survey 
is deep enough that the bulk of the stellar Ught over this red- 
shift range is directly observed. As far as possible, we comple- 
ment our data by "local" information from low-redshift surveys 
(e.g. 2dFGRS and SDSS), as our survey volume is too small for 
z < 0.2. 

The observed galaxy SEDs cover a wavelength range of 
Aobs = 350 . . . 930 nm and allows us to determine restframe 
properties of the population at 280 nm and at B with hardly 
any extrapolation in the range 0.2 < z < 1.2. We also derived 
restframe properties in the SDSS r-band, by applying a type- 
dependent extrapolation of the observed SED at z ^ 0.5. 

The goal of this present paper is to present an empirical pic- 
ture of the evolution of the galaxy population as a function of 
luminosity, SED and redshift, as a solid constraint for any mod- 
els of statistical galaxy evolution. Direct comparison with cos- 
molgical models will be reserved for future papers. Our main 
findings are as foUows: 

1. Survey areas approaching an area of 1 are necessary to 
reduce field-to-field variations at L* below a factor of two 
in redshift bins of Az ^ 0.2. 

2. The faint end slope of the luminosity function depends 
quite strongly on SED-type at all observed redshifts. 
However, within a given (non-evolving) SED type there is 
no evidence for evolution of the luminosity function shape 
with redshift (to z 1). 

3. A comparison of our "quasi-local" sample at z = [0.2, 0.4] 
with 2dFGRS and SDSS shows consistent results around 
the knee of the luminosity function, but a steeper faint end 
in COMBO-17, probably reflecting the evolution of faint 
blue galaxies between the median redshifts ^ 0.05 and ~ 
0.34. 

4. If we define SED type via a fixed, present-day spectral se- 
quence for galaxies, from ellipticals to young starbursts, 
then the evolution of the luminosity function from z ~ 1.2 
to now depends strongly on SED type, both in its density 
normalization, (j)*, and in its characteristic luminosity, M*: 

(a) The faint end slope isa = 0.50±0.20for early type 
galaxies (SED typel), steepening to a = 1.50 ± 0.06 for 
galaxies with the bluest starburst colours (SED typeA-). 

(b) For early type galaxies (defined by z=0 colours bluer 
than Sa galaxies) </>* increases by more than an order of 



magnitude from z = 1.2 to now, while M* gets fainter by 
more than a magnitude. 

(c) For the latest type galaxies ((>* remains roughly constant 
over the redshift range probed, while gets fainter by 
almost two magnitudes. 

5. When averaging over all SED types, the evolution of 
the luminosity density depends considerably on the rest- 
wavelength considered: 

(a) In the restframe B- and r-bands, the integrated luminos- 
ity density remains constant from z = 1.1 to z = 0.5, drop- 
ping subsequently to the present epoch by perhaps 30%. 

(b) In the restframe near-UV (280 nm) the integrated lu- 
minosity density drops by a factor of six from z ~ 1 to 
now, where it appears that much of that drop occurs at red- 
shifts below z ^ 0.6. Our more accurate estimates of the 
near-UV luminosity density imply a considerably shallower 
evolution than indicated by the data of Lilly et al. (1995), 
but is consistent within their 1. 5a confidence limits. 

6. The ratio of the stellar luminosity density contributed by 
early and late SED type galaxies at different wavelengths 
changes quite drastically with redshift: At redshift z 1 
strongly star-forming galaxies itype3 and 4) make up the 
bulk of the radiation in aU three restframe bands with a 
contribution near 90% at 280 nm. At the present epoch, in 
particular the extreme starbursts (typeA) play only a signifi- 
cant role at 280 nm, providing about half of the luminosity 
there. At the same time, the fraction of light from galaxies 
with early- type spectra (typel) have increased from below 
10% to 40% in the B- and r-bands over the same period. 

We reserve any rigorous modehng of these results for a later 
paper, as it appears that the data warrant a cosmologically mo- 
tivated interpretation, such as a comparison with semi-analytic 
or N-Body/SPH models rather than any empirical speculation 
about "luminosity evolution" or "density evolution" of various 
sub-populations. Nonetheless, it seems approriate to make a 
few quaUtative remarks in closing: our results confirm earUer 
findings that the epoch of galaxy formation, or star formation 
in galaxies, has been waning over the second half of the present 
cosmic age, leading to a strong decrease of star-burst galax- 
ies (both luminous and faint) accompanied by the strong emer- 
gence of "old" galaxies, whose overaU SED is redder than an Sa 
template spectrum. The COMBO-17 data set permits to study 
this evolution differentially for different SED types and still 
provides much smaller error bars, due to the enlarged survey 
area and sample increase. Compared to earlier work, the de- 
cline of cosmic star formation activity appears gentler, i.e. the 
near-UV luminosity density as a proxy for the density of young 
stars, is found to drop only by a factor of five from z = 1.1 to 
now. 

Beyond questions of interpretation, there are a number of 
observational issues that still need to be addressed. Calibration 
of our 17-band redshifts through spectra for samples of many 
hundred objects is forthcoming in the course of the ESO- 
GOODS key project, and will increase the reUability of our 
estimates even further. The addition of several other fields wiU 
reduce the impact of the field-to-field variation, now still a hm- 
iting factor in several aspects of the analysis. Finally, the full 
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analysis of the SEDs and the addition of near-IR data will al- 
low us to extend the analysis to include direct observational 
estimates of the stellar mass. 
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Appendix A: Luminosity functions piots and 
data 

The appendix contains detailed tables and plots of luminos- 
ity functions for all restframe passbands, all SED types and 
all redshift intervals. Whenever, the reader likes to refer to a 
parametric luminosity function for all galaxy types combined 
we strongly suggest to use a sum of the type-specific functions, 
rather than the Schechter fit to the combined sample, which is 
not entirely appropriate, especially in the extrapolated faint do- 
main. For this reason, the luminosity density of all galaxy types 
combined is given as the sum of the four types rather than cal- 
culated from the inappropriate Schechter fit. 

Table A.l. Field-to-field variation of <}>* in the "quasi-local" 
sample of 5,674 galaxies at ^; = [0.2, 0.4]. 
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53.27±12.98 


280 


100.02 
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184.75 


158.41±50.65 



28 



Wolf et al.: The galaxy luminosity function from COMBO-17 



LF redshift evolution in r-band for SED TYPE 1 with a =+0.43 
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Fig. A.l. Redshift evolution of 0(A/r) for type 1 galaxies: VJnax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black fine, while the fits for the other redshifts are shown as grey lines. Local 
reference: SDSS data. 

Table A.2. Redshift evolution of parameters for (j){M^): Results of STY fits in five redshift intervals centered at z — 
[1.1,0.9,0.7,0.5,0.3]. The faint-end slope a was determined in the quasi-local sample and fixed for all redshifts. Listed are 
M*, (j)* with its field-to-field variation, j and the covariance between 0* and L* . The j values for the combined sample are 
actually sums over the jtype- 



parameter 


< z> 


M' - 5 log h 
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(Vega mag) 
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0.3 


-20.92±0.14 


44.09±22.81 


0.43±0.18 


8.04±4.16 


-0.009 




0.5 
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Fig. A.2. Redshift evolution of 0(A/r) for type 2 galaxies: V;„ax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black fine, while the fits for the other redshifts are shown as grey lines. Local 
reference: SDSS. 
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LF redshift evolution in r-band for SED TYPE 3 with a = -1.33 




Fig. A.3. Redshift evolution of 0(A/r) for type 3 galaxies: VJnax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black fine, while the fits for the other redshifts are shown as grey lines. Local 
reference: SDSS. 
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LF redshift evolution in r-band for SED TYPE 4 with a = -1.50 




Fig. A.4. Redshift evolution of 0(A/r) for type 4 galaxies: VJnax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black fine, while the fits for the other redshifts are shown as grey lines. Local 
reference: SDSS. 



32 



Wolf et al.: The galaxy luminosity function from COMBO-17 



LF redshift evolution in r-band for ALL SED TYPES with a=free 
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Fig. A.5. Redshift evolution of 4>{M,-) for all galaxies: Vmax data points are shown with error bars for one redshift interval per 
panel. The corresponding STY fit is plotted as a black line, while the fits for the other redshifts are shown as grey lines. Local 
reference: SDSS. 
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LF redshift evolution in B-band for SED TYPE 1 with a =+0.52 
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Fig. A.6. Redshift evolution of <f>{MB) for type 1 galaxies: Vmax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black line, while the fits for the other redshifts are shown as grey lines. Local 
reference: 2dFGRS. 



Table A.3. Redshift evolution of parameters for (/)(Mb): Results of STY fits in five redshift intervals centered at z = 
[1.1,0.9,0.7,0.5,0.3]. The faint-end slope a was determined in the quasi-local sample and fixed for all redshifts. Listed are 
M*, </)* with its field-to-field variation, j and the covariance between (/>* and L*. The j values for the combined sample are 
actually sums over the jtype- 
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< z > 


M' - 5 log /i 
(Vega mag) 


4>* X lo^'* 


a 
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(h/Mpc^) 






0.3 


-19.49±0.15 


48.68±15.96 


0.52±0.20 


5.41±1.77 


-0.049 




0.5 


-19.79±0.16 


28.49±12.04 


± 


4.17±1.76 


-0.035 


type 1 


0.7 


-20.19±0.16 


13.44± 3.62 


± 


2.85±0.77 


-0.144 




0.9 


-20.06±0.21 


4.63± 4.16 


± 


0.87±0.78 


-0.788 




1.1 


-20.34±0.23 


2.46± 1.99 


± 


0.60±0.48 


-1.363 



0.3 


-19.72±0.17 


40.97±19.05 


-0.90±0.08 


3.97±1.85 


-0.468 


0.5 


-20.07±0.19 


49.56± 6.48 


± 


6.63±0.87 


-0.407 


0.7 


-20.25±0.11 


48.67± 2.43 


± 


7.72±0.39 


-0.646 


0.9 


-20.44±0.12 


38.48±20.29 


± 


7.24±3.82 


-1.151 


1.1 


-20.75±0.15 


18.63± 6.00 


± 


4.66±1.50 


-1.522 



0.3 


-20.11±0.23 


19.30± 


9.45 


-1.34±0.07 


3.84±1.88 


-0.684 


0.5 


-20.52±0.37 


28.44± 


7.33 


± 


8.25±2.13 


-0.625 


0.7 


-20.38±0.16 


27.61± 


4.31 


± 


7.04±1.10 


-0.836 


0.9 


-20.50±0.14 


31.43± 


8.82 


± 


8.95±2.51 


-1.149 


1.1 


-20.60±0.15 


23.51± 


2.94 


± 


7.34±0.92 


-1.371 



0.3 


-19.05±0.16 


63.44± 4.93 


-1.47±0.06 


5.76±0.45 


-1.030 


0.5 


-19.45±0.17 


69.93±14.50 


± 


9.24±1.92 


-0.914 


0.7 


-19.44±0.12 


62.25±12.51 


± 


8.11±1.63 


-1.278 


0.9 


-19.65±0.14 


64.13±14.68 


± 


10.15±2.32 


-1.614 


1.1 


-20.08±0.12 


39.11±13.04 


± 


9.17±3.06 


-1.672 





0.3 


-20.50±0.13 


53.27±12.98 


-1.47±0.03 


18.97±3.21 


-0.721 




0.5 


-20.03±0.06 


179.34±31.08 


-1.09±0.03 


28.30±3.47 


-0.526 


all 


0.7 


-21.00±0.11 


41.74± 3.88 


-1.63±0.04 


25.72±2.15 


-1.028 




0.9 


-21.27±0.09 


23.34±10.91 


-2.00±0.02 


27.21±5.18 


-1.672 




1.1 


-21.20±0.09 


33.81±10.84 


-2.00±0.03 


21.78±3.56 


-2.225 




Fig. A.7. Redshift evolution of <f>{MB) for type 2 galaxies: Vmax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black line, while the fits for the other redshifts are shown as grey lines. Local 
reference: 2dFGRS. 
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LF redshift evolution in B-band for SED TYPE 3 with a = -1.34 




Fig. A.8. Redshift evolution of <f>{MB) for type 3 galaxies: Vmax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black line, while the fits for the other redshifts are shown as grey lines. Local 
reference: 2dFGRS. 




Fig. A.9. Redshift evolution of <f>{MB) for type 4 galaxies: Vmax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black line, while the fits for the other redshifts are shown as grey lines. Local 
reference: 2dFGRS. 
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Fig. A. 10. Redshift evolution of 0(A/b) for all galaxies: Vmax data points are shown with error bars for one redshift interval per 
panel. The corresponding STY fit is plotted as a black line, while the fits for the other redshifts are shown as grey lines. Local 
reference: 2dFGRS. 



38 



Wolf et al.: The galaxy luminosity function from COMBO-17 




Fig. A.ll. Redshift evolution of (j){M2so) for type 1 galaxies: Knax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black Une, while the fits for the other redshifts are shown as grey lines. No 
local reference available. 

Table A.4. Redshift evolution of parameters for 0(Af28o): Results of STY fits in five redshift intervals centered at z = 
[1.1,0.9,0.7,0.5,0.3]. The faint-end slope a was determined in the quasi-local sample and fixed for all redshifts. Listed are 
M*, </)* with its field-to-field variation, j and the covariance between </)* and L*. The j values for the combined sample are 
actually sums over the jtype- 



parameter 


< z > 


M* - 5 log h 


<p' X 10"" 


a 


j X lO^Lo 








(Vega mag) 


(h/Mpcr^ 




(h/Mpc^) 






0.3 


-17.81±0.14 


39.49±19.00 


0.56±0.19 


3.36± 1.62 


-0.413 




0.5 


-18.49±0.16 


23.19± 9.20 


± 


3.69± 1.47 


-0.176 


type 1 


0.7 


-18.49±0.14 


14.90± 5.15 


± 


2.38± 0.82 


-0.175 




0.9 


-18.76±0.19 


4.09± 2.54 


± 


0.84± 0.52 


-0.599 




1.1 


-19.34±0.23 


3.09± 2.05 


± 


1.08± 0.72 


-1.321 




0.3 


-19.58±0.20 


20.84±11.66 


-1.20±0.05 


7.63± 4.27 


-0.687 




0.5 


-20.03±0.29 


29.04± 6.57 


± 


16.09± 3.64 


-0.599 


type 2 


0.7 


-19.70±0.11 


34.23± 2.25 


± 


14.00± 0.92 


-0.862 




0.9 


-20.08±0.13 


24.12± 7.81 


± 


13.94± 4.51 


-1.130 




1.1 


-20.51±0.17 


16.73± 1.11 


± 


14.43± 0.95 


-1.510 




0.3 


-20.43±0.25 


11.75± 5.51 


-1.48±0.04 


13.64± 6.39 


-0.741 




0.5 


-21.07±0.44 


17.33± 4.04 


± 


36.15± 8.44 


-0.695 


type 3 


0.7 


-20.60±0.16 


24.44± 3.24 


± 


33.19± 4.39 


-0.949 


0.9 


-20.83±0.13 


23.29± 7.79 


± 


39.29±13.14 


-1.216 




1.1 


-20.96±0.14 


20.72± 4.18 


± 


39.23± 7.92 


-1.412 




0.3 


-19.39±0.14 


74.43±14.93 


-1.53±0.05 


36.17± 7.26 


-0.977 




0.5 


-20.09±0.20 


65.77±14.89 


± 


61.21±13.86 


-0.990 


type 4 


0.7 


-20.01±0.11 


68.99±13.12 


± 


59.90±11.39 


-1.419 




0.9 


-20.20±0.09 


89.55±24.69 


± 


92.73±25.57 


-1.823 




1.1 


-20.81±0.11 


48.60±23.77 


± 


88.54±43.31 


-1.729 




0.3 


-19.41±0.08 


158.41±50.65 


-1.30±0.03 


60.80±10.69 


-0.773 




0.5 


-19.41±0.Q5 


374.96±71.17 


-0.72±0.04 


117.16±16.69 


-0.656 


aU 


0.7 


-19.80±0.06 


193.29±29.69 


-1.10±0.05 


109.47±12.27 


-1.221 




0.9 


-20.31±0.10 


148.33±45.46 


-1.34±0.09 


146.80±29.10 


-1.580 




1.1 


-20.80±0.16 


88.56±28.57 


-1.47±0.17 


143.29±44.05 


-1.694 
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Fig. A.12. Redshift evolution of (j){M2so ) for type 2 galaxies: Knax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black hne, while the fits for the other redshifts are shown as grey lines. No 
local reference available. 
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LF redshift evolution in 280-band for SED TYPE 3 with a = -1.48 




Fig. A.13. Redshift evolution of (j){M2so ) for type 3 galaxies: Knax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black hne, while the fits for the other redshifts are shown as grey lines. No 
local reference available. 
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LF redshift evolution in 280-band for SED TYPE 4 with a = -1.53 




M280 - 5 log h Maao - 5 log h M280 - 5 log h 



Fig. A.14. Redshift evolution of (j){M2so ) for type 4 galaxies: Knax data points are shown with error bars for one redshift interval 
per panel. The corresponding STY fit is plotted as a black hne, while the fits for the other redshifts are shown as grey lines. No 
local reference available. 
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LF redshift evolution in 280-band for ALL SED TYPES with a=free 




Fig. A.15. Redshift evolution of 0(M28o) for all galaxies: Vinax data points are shown with error bars for one redshift interval per 
panel. The corresponding STY fit is plotted as a black line, while the fits for the other redshifts are shown as grey lines. No local 
reference available. 
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